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(57) Abstract: A variable density non-Cartesian parallel-imaging method for reconstructing a magnetic resonance (MR) image is 
provided. In embodiments of the invention, an MR data set is obtained by sampling first and second sampling regions, wherein a 
first region is sampled with a first sampling density that is higher than a second sampling density of a second region. MR images 
corrupted by aliasing artifacts are reconstructed from the data obtained with each one of the coil-elements of a coil array. These 
images can be combined into one, de-aliased image using a modified version of Cartesian SENSE. The modification allows all the 
available k-space lines to be used in the processing, despite the fact that different k-space regions have different sampling densities 
(i.e. non-Caitesian sampling). Using all available lines is advantageous in terms of signal-to-noise ratio. Advantages of embodiments 
of the invention over previous methods also able to deal with non-Cartesian sampling schemes may include one or more of simplicity, 
ease of implementation, not having to fit sensitivities to taiget functions as part of the reconstruction, fast processing speed and/or 
the avoidance of possible eirors resulting from solving laige systems of equations. 
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VARXABLE K-SPACE-DENSITY PARALLEL MAGNETIC RESONANCE IMAGING 

Field of the Invention 
The invention relates generally to reconstruction methods for magnetic 
5 resonance unage (MRI) data and more specifically to reconstruction methods for MRI 
data acquired using non-Cartesian sampling and parallel imaging. 

Background of the Invention 
In many magnetic resonance imaging (MRI) applications, a given region of 

10 tibe body is imaged repeatedly to capture its time variations. For example, such 

dynamic applications include functional MRI (in which brain changes are mduced by 
a time-varying paradigm), time-resolved angiography (in which changes in the blood 
vessels are caused by tiie passage of a bolus of contiast agent), and cardiac imaging 
(in vMch the heart changes as it beats, and also possibly as a bolus of contrast ^ent 

1 5 passes through it). The temporal resolution of tiie MRI, i.e. tiie time to acquire a time 
fi-ame, should be good enoi^ to capture the important features of the temporal 
changes. In tiie event that tiie readily available temporal resolution proves 
insufficient, tiiere exist many dynamic MRI mefliods able to improve it. Some of 
tiiese metiiods include UNFOLD, parallel imaging (e.g. SMASH, SENSE, SPACE- 

20 RIP) and partial-Fourier imaging techniques. Through some assumption(s) and/or tiie 
use of prior information, tiiese metiiods allow a fraction of tiie required data to be 
calculated instead of measured. This reduction in tiie amount of acquired data usually 
translates diiectiy into a corresponding reduction in the time to acquire the data and 
thus can improve the temporal resolution, and/or the total scan time. 

1 
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Possibly due to its simplicity and fast processing speed, Cartesian SENSE may 
be the most commonly used parallel imaging method. Cartesian SENSE is typically 
limited to Cartesian, regular sampling schemes. More general methods such as 
SPACERIP or general SENSE can be used to reconstruct data acquired along more 

5 complicated trajectories in k-space. Recently, non-Cartesian sampling schemes such 
as variable-density SMASH, GRAPPA and others have allowed sensitivity 
information to be obtained as part of the dynamic acquisition by sampling more 
densely the center of k-space than the outer regions. These "self-calibrated" methods 
do not require the acqmsition of a reference scan to measure the coil sensitivity; 

1 0 instead, coil sensitivity is preferably calculated directly from the fully sampled region 
around the center of k-space, in the dynamically acqxiired data set 

A regular Cartesian sampling scheme 10 is shown in Fig. lA in which only one 
k-space line 12 out of every four is sampled. Fig. IB shows the resulting image 20 
from one of the coils. The image 20 is corrupted by aliasing artifacts resulting from 

1 5 sampling only 24 k-space lines rather than 96. As described in the publication 
SENSE: sensitivity encodmg for fast MRI, 42 Magn Reson Med 952 (1999) by 
Klaas P. Pruessmann et al, Cartesian SENSE can be used to separate the four 
overlapped spatial locations at each image pixel to reconstruct a fidl, de-aliased 
image. Fig. IC shows an image 30 that was reconstructed and de-aliased using 

20 Cartesian SENSE. 

Fig. 2A shows another sampling technique 40 that is similar to the technique 
10 used in Fig. 1 A with the exception that all of the missing k-space locations 42 are 
mied with zeros before applying a Fourier transform (FT). As a consequence of the 



2 
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presence of the zeros, a full field of view (FOV) is reconstructed in the image 50 
shown in Fig. 2B instead of the smaller, acquired FOV of the image 20 m Fig. IB. 

In an alternative method, the four overlapped spatial locations at each pixel of 
the image 50 in Fig. 2B covdd be separated by keepmg only the non-aliased pixel and 

5 discarding the remaining three aliased pixels. This method and the method described 
above for Fig. IB are mathematically equivalent, and lead to numerically identical 
treated images 30 and 60 as shown in Figs. IC and 2C. 

Typically, it is deskable to use a Cartesian sampling function because of the 
simplicity of calculations associated with the data reconstruction. Other known non- 

10 Cartesian sampling methods, such as SPACERff and the general version of SENSE 
allow data acquired with a non-Cartesian sampling scheme to be reconstructed into 
images, but require a significantly larger number of calculations to do so. While 
variable-density SMASH and GRAPPA do allow the use of a sampling strategy that 
departs fi:om a Cartesian grid, these methods make the approximation that coil 

1 5 sensitivities can be combined to emulate functions that are related to Fourier basis 
functions, an approximation absent in the SENSE/SPACERIP approach. 

Summary of the Invention 

A variable density, non-Cartesian method for reconstructing a magnetic 
20 resonance (MR) image is provided. In an embodiment of the invention, an MR data 
set is obtained by sampling first and second sampling regions using parallel imaging, 
wherein a first region is sampled with a first sampling density that is higher than a 
second sampling density of a second region. The MR data set can be reconstructed 
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using a modified method of Cartesian SENSE that includes reconstructing both the 
first and the second sampling regions using an acceleration rate of the second region. 

In embodiments of the invention, a method of magnetic resonance (MR) 
imaging for leconstructing an MR image is provided, wherein the method generally 
5 includes: acquiring an MR data set by using a phased-array coil consisting of a 
plurality of individual coil-elements and by sampling first and second sampling 
regions, wherein the first region is sampled with a higher sampling density than the 
second region; reconstructing a plurality of images fi:om a plurality of individual coil- 
elements ; and combining the plurality of unages from the individual coil-elements 
10 mto a resultant MR image by using parallel-imagmg reconstruction. 

In embodiments of the invention, the sampling density is higher near k-space 
center than elsewhere, and the data is reconstructed using a variable-density version 
of Cartesian SENSE. 

Embodiments of the invention reconstruct MR images using a number of 
15 calculations almost as low as that of Cartesian SENSE, while allowing k-space 
sampling schemes that depart from purely Cartesian grids. 

Various embodiments of the invention may provide one or more of the 
following advantages. Advantages of embodiments of the invention over previous 
methods also able to deal with non-Cartesian sampling schemes may include one or 
20 more of: simplicity in data reconstruction, ease of implementation (similar to 
Cartesian SENSE), not having to fit sensitivities to target functions as part of the 
reconstruction, fast processing speed and the avoidance of possible errors resulting 
from solving large systems of equations. Advantages of the invention may include a 
faster processing speed than non-Cartesian methods such as SPACERIP and general 

4 
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SENSE, and the absence of the need for fitting sensitivity maps to target functions as 
in variable-density SMASH and GRAPPA. Non-Cartesian sampling can be used and 
im^es reconstructed without using large numbers of calculations. 



5 
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Brief Description of the Figures 
Fia 1 A is a diagram of a prior art Cartesian SENSE sampling scheme having 
acceleration rate of four. 

FIGS. IB and IC are MR images corresponding to the data acquired in FIG. 
5 1 A, wherein FIG. 1 B is an image from one of a plurality of individual coil-elements 
of a phased-array coil, and FIG. IC is the resultant image after Cartesian SENSE 
processing. 

FIG. 2 A is a diagram of a prior art Cartesian SENSE sampling scheme similar 
to the sampling scheme of FIG. 1 A but with zeros filled in missing k-space locations. 
10 FIGS. 2B and 2C are MR images corresponding to the data acquired in FIG. 

2A, wherein FIG. 2B is an image fi:om one of a plurality of individual coil-elements 
of a phased-array coil and FIG. 2C is the resultant image after Cartesian SENSE 
processing. 

FIG. 3 A is a diagram of a sampling scheme according to the invention. 
15 FIGS. 3B and 3C are MR images corresponding to the data acquired in FIG. 

3 A, wherein FIG. 3B is an image &om one of a plurality of individual coil-elements 
of a phased-array coil and FIG. 3C is the result image fi:om an embodiment of the 
invention. 

FIG. 4 is a simplified diagram of a phased array of coils for MR imaging an 

20 object. 

FIG. 5 is a simplified schematic diagram of a system for unaging an object in 
accordance witih the invention. 

FIG. 6 is a block flow diagram of a process of imaging an object usmg the 
system shown in FIG. 5. 
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Detailed Description of the Invention 
Embodiments of the mvention provide a method for modifying the MR 
reconstruction technique of Cartesian SENSE such that it can acconmiodate a 

5 sampling strategy that departs j&om a purely rectangular grid and allows for MR data 
reconstruction with a relatively low number of calculations. Near-Cartesian data 
sampling can be performed in conjunction with parallel processing of acquired data. 
The near-Cartesian sampling departs from pure Cartesian sampling, although more 
general methods such as SP ACERIP and GRAPPA can handle more drastic 

10 departures from Cartesian sampling. Data are more densely sampled near a center of 
k-space than away from the k-space center, and receiver-coil sensitivity is measured 
"on the fly" during data acquisition to provide a self-calibratmg system. 

For purpose of illustration, the invention will now be described with reference 
to Figs. 3 A-3C, in which a method for reconstructing MR image data is shown. Fig. 

1 5 3 A shows an exemplary sampling method 70 in which a first region 72 of k-space is 
sampled at a higher density than the remaining sampled regions 74 and 76. Region 
72, which includes the center of k-space, is sampled at an acceleration rate of 2 such 
that one out of every two k-space Unes 77 is sampled. Regions 74 and 76, which are 
at the upper and lower outer edge regions of k-space, respectively, ate sampled at an 

20 acceleration rate of 4 such that one out of every four k-space lines 78 is sampled. 

The image shown in Fig 3B, which has desired information 91 and aliasing 
artifacts 92, 94 and 96, is obtained by correctmg for sampling density and Fourier 
transforming the data obtained from one of the individual coil-elements fonning the 
phased-array coil while using the sampUng strategy of Fig. 3A. In Fig. 3B, aliasmg 

7 
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artifacts 92, 94 and 96 are satisfectory for use with embodiments of the invention as 
they look fairly sharp and are not diffused or smeared over the image. Such non- 
smeared artifacts can be expected if the non-Cartesian sampling strategy mildly 
departs from a pxirely Cartesian strategy (as is the case in FIG. 3 A). To obtain the 

5 reconstructed image 100 shown in Fig. 3C, the method of Cartesian SENSE modified 
according to the invention is applied to remove the aliasing artifacts 92, 94 and 96. 
That embodiments of the invention can remove the aliasing artifacts 92, 94 and 96 
similar to Fig. IC and 2C indicates that the invention can be used to modify Cartesian 
SENSE to handle small departures from regular, Cartesian sampling. For example, 

10 the invention can be used to increase the sampling density near k-space center, as 
proposed in other methods for self-navigated parallel imaging. 

As part of the invention, a region of k-space is sampled to obtain data at a 
higher k-space sampling density and a lower acceleration rate than another region and 
the sampled data is reconstructed using a modified version of Cartesian SENSE. For 

15 example, a method of the invention could include sampling at a higher density near a 
k-space center than elsewhere, allowing self-calibration and/or a reduction in artifact 
content. Alternatively, it may be more desirable to sample a region other than the 
center of k-space at a higher sampling rate. The invention allows all of the lines 
acquired in FIG. 3 A (i.e. a non-Cartesian set) to contribute to the reconstruction, 

20 which is desirable in terms of signal-to-noise-ratio (SNR). One could dismiss lines in 
FIG. 3 A to obtain a Cartesian pattern, and reconstruct the resulting Cartesian data set 
with Cartesian SENSE. However, dismissing lines comes at a cost in SNR, which can 
be avoided with the invention. In the reconstruction, the acceleration of the least 
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densely sampled region, i.e. the highest acceleration of all the sampled k-space 
regions, is used in the parallel imaging algorithm. 

Typically, as part of the method of reconstructing the MR data, a system of 
linear equations will be solved during which two steps will be performed. These steps 
5 1) invert a matrix, and 2) multiply the inverted matrix by a particular vector. Step 1 
uses most of the processing time and step 2 is performed much quicker. Step 1 is 
represented in the above-mentioned article by Pniessmann et al by approximating the 
inverse of a sensitivity matrix S given by: 

10 where Y,p counts coils (see FIG. 4 and description below) and the superimposed 
pixels, respectively, rp denotes the position of the pixel p, and Sy is the spatial 
sensitivity of the coil y. The inverse of S is an unfolding matrix U given by: 

U = (SV^S)-^SV^ 
where the superscript H indicates the transposed complex conjugate, and ^ is a 
15 receiver noise matrix. In the second step, the unfolding matrix U is multiplied by a 
vector a, that represents the complex image values of a chosen pixel, to obtain a 
resulting vector v given by: 
v = Ua. 

For step 1, an exemplary embodiment of the invention reconstructs the data 
20 using the higher acceleration factor of two regions of k-space having different 

acceleration factors. Suppose, for example, an image is generated which consists of N 
pixels, in which a phased-array coil made of a collection of M different individual 
coils is used, and that parallel imaging is used to speed up data acquisition by a factor 
of n. The mefliod of Cartesian SENSE uses the inversion of N/n matrices with each 
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matrix having up to n*M elements. Embodiments of the invention slightly differ from 
Cartesian SENSE in that the higher acceleration factor (nh) of the two regions is used 
for calculating the nimiber of matrices (N/nn) and the number of elements of each 
matrix (up to nh*M). The total acceleration factor n, v^rhich takes the whole sampling 
5 scheme into account (and not just the region sampled with acceleration Uh), is smaller 
thanuh. 

In step 2, an exemplary embodiment of the invention mxiltiplies the inverted 
matrices "Uh" more times than is done in Cartesian SENSE (N multiplications rather 
than N/Uh, i.e. one multiplication per pixel in the reconstructed field-of-view rather 

1 0 than one multiphcation per pixel m the smaller, acquired FOV as required in 
Cartesian SENSE). Because this step 2 tends to be much faster than step 1, 
embodiments of the invention thus use marginally more processing than Cartesian 
SENSE, and very significantiy less than other variable-density methods such as 
SPACERIP and general SENSE (as those use the inversion of much larger matrices 

1 5 than in Cartesian SENSE and in the invention). 

Embodiments of the invention are used to handle only relatively small 
departures from a Cartesian, regular samplmg. In exemplary preferred embodiments, 
a regular Cartesian sampling grid in k-space, e.g., the technique 10 shown in FIG. 1 A, 
is perturbed using another regular or substantially regular pattern. Preferably, the 

20 sampling in k-space is not modified by a random insertion of k-space lines. For 
example, the sampling scheme depicted m FIG. 3A consists of two uiterleaved 
Cartesian sampling schemes: one line every four over all of k-space, and one Ime 
every four over the center of k-space and shifted by one line with respect to the 
former. Non-regular sampling schemes might lead to smeared artifacts, for which 

10 
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SPACERIP or general SENSE could be used. Using the invention, the near-Cartesian 
sampling will produce well-localized replica of the object which can be removed . 
In exemplary embodiments of the invention, missing k-space locations 80 from the 
sampling method 70 shown in Fig. 3A are filled with zeros and corrected for sampling 

5 density. With these embodiments a conventional Fast-Fourier transform algorithm 
may be used to generate the images from individual coil-elements. 

Refening to FIG. 4, an arrangement for MR imaging includes an array 1 10 of 
coils 1 12. The coils 1 12 are configured to acquire data of an object 1 14, here a 
person's torso including a heart 1 16. Measurements obtained by the coils 1 12 are 

1 0 relayed through cables 1 1 8 to appropriate processing devices. 

Referring to FIG. 5, an imaging system 120 includes the coils 1 12, receiver 
boards 122, and a processor 124. An MR scanner can be used to generate MR images 
for the invention. One such MR scanner is manufactured by General Electric under 
the brand name Signa. The MR scanner can include the phased array 1 1 0 of coils 

1 5 having a plurality of different coil-elements, or the phased array 1 1 0 of coils can be 
separate firom the MR scanner. The receiver boards 122 receive the measured data 
from the coils 1 12 and transmit these data to tiie processor 124. The processor 124, 
such as that found in a standard computer or a laptop, can be used to combine 
generated MR images. Alternatively, tiie processor 124 can be included as part of the 

20 MR scanner. 

Methods of the invention can be performed with a computer and computer 
software installed thereon that has instructions to perform steps of the invention. 
Alternatively, methods of the invention can be performed with equipment that has 
installed hardware or firmware having instructions to perform steps of the invention, 

11 
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Software used with embodiments pf the invention can be stored on computer usable 
media for storing data, such as, for example, but not limited to, floppy disks, magnetic 
tape, zip disks, hard drives, CD-ROM, ROM, optical disks, or a combination of any of 
these. 

5 In operation, referring to FIG. 6, with further reference to FIGS. 4-5, a process 

130 of imaging the object 1 14 using the system 120 includes the stages shown. The 
process 130, however, is exemplary only and not limiting. The process 130 can be 
altered, e.g., by having stages added, removed, or rearranged. 

At stage 132, raw image data are obtained. The coils 112 measure data from 
1 0 the object 1 1 4, e.g., including the heart 1 1 6. These data are relayed through the 
cables 1 1 8 to the receiver boards 122 and on to the processor 124. The data are 
obtained by the coils 1 12 using a near-Cartesian sampling scheme, such as that shown 
in FIG. 3 A. 

At stage 134, the processor 124 produces an image of the object 1 14 from 
15 each one of the receiver coils, the image including aliasing effects. The processor 124 
uses data from the coils near the center of k-space to self-calibrate the system 120 by 
determining the sensitivities of the coils 1 12. 

At stage 136, using the images from the individual receiver-coils and the 
determmed sensitivities, the processor 124 removes the aliasing artifacts from the 
20 image. The processor 124 uses parallel processing techniques and knowledge of the 
receiver-coil sensitivity to remove the aliasing artifacts and yield an image 
substantially free of the artifacts while using a relatively low amount of calculations. 

Advantages of embodiments of the invention over previous methods also able 
to deal with non-Cartesian sampling schemes include one or more of simplicity, ease 

12 
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of implanentation, not having to fit sensitivities to taj^et fiinctions as part of the 
reconstruction, fast processing speed and the avoidance of possible errors resulting 
from solving large systems of equations. 

Having thus described at least one illustrative embodiment of the invention, 
various alterations, modifications and improvements will readily occur to those 
skilled in the art. Such alterations, modifications and improvements are intended to 
be within the scope and spirit of the invention. Accordmgly, the foregoing 
desCTiption is by way of example only and is not mtended as limiting. The 
invention's limit is defined only in the following claims and the equivalents thereto. 
What is claimed is: 



13 
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1 . A method for reconstructing a magnetic resonance (MR) image, the method 
comprising: 

obtaining an MR data set by sampling first and second sampling regions, 
wherein a first region is sampled with a first sampling density that is higher than a 
5 second sampling density of a second region; 

generating MR images corrupted by aliasing artifacts by reconstructing the 
obtamed MR data set; and 

combining the generated MR images with parallel imagmg, using an 
acceleration rate of the second region. 

10 

2. The method of claim 1 , wherein obtaining an MR data set includes acquiring 
the MR data set using a phased-array coil consisting of a plurality of coil-elements, 
and wherem combining the generated MR images includes combining a plurality of 
images bom a plurality of individual coil-elements into a resultant image. 

15 

3. The method of claim 2, wherein combining the plurality of images firom the 
individual coil-elements includes invertmg a number of matrices, the number of 
matrices being equal to a number of image pixels in the resultant MR image divided 
by the acceleration rate of the second region, wherein each matrix has a number of 

20 elements equal to or less than the number of coil-elements in the phased array coil 
times the acceleration rate of the second region. 

4. The method of claim 2, wherein combining a plurality of images includes 
filling missing k-space locations with zeros. 

14 
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5. The method of claim 4, further comprising the application of a Fourier- 
transform after filling missing k-space locations with zeros. 

5 6. The method of claim 1, wherein the first sampling region is substantially at a 
center of k-space. 

7. The method of claim 1 wherein the first and second sampling densities both 
provide substantially regular sampling patterns. 

10 

8. The method of claim 7 wherein the first samphng density provides a first, 
Cartesian sampling pattern and the second sampling density provides a second 
sampling pattem comprising a Cartesian pattern perturbed by a substantially regular 
pattern. 

15 

9. The method of claim 1 wherein the sampling is performed in a near-Cartesian 
manner. 

10. A method of magnetic resonance (MR) imaging for reconstructing an MR 
20 image, the method comprising: 

acquiring an MR data set by using a phased-array coil consistmg of a plurality 
of individual coil-elements and by sampling first and second sampling regions, 
wherein the first region is sampled with a higher sampling density than the second 
region; 

15 
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reconstructing a plurality of images from a plurality of individual coil- 
elements; and 

combining the plurality of images from the individual coil-elements into a 
resultant MR unage using parallel imaging. 

5 

1 1 . The method of claim 1 0, \^erein combining the plurality of images includes 
inverting a number of matrices, the number of matrices being equal to a number of 
image pixels in the resultant MR image divided by the accelerated rate of the second 
region, wherein each matrix has a number of elements equal to or less than the 

1 0 number of independent coil-elements of the phased array coil times an accelerated 
rate for the second region. 

12. The method of claim 10, wherein combining the plurality of images includes 
inverting more than N matrices, each having less than NNM/n elements, wherem N is 

1 5 equal to a number of image pixels in one or more phase-encoded directions for a 
location in a frequency-encoded direction in the resultant MR image, M is equal to a 
number of individual coil-elements in the phased array coil and n is equal to an 
accelerated rate for the second region. 

20 13. The method of claim 1 0, wherem reconstructing a plurality of images from the 
individual coil-elements includes filling missing k-space locations with zeros and . 
applying a Fourier-transform. 



16 
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14. The method of claim 13, wherein the first sampling region is substantially at a 
center of k-space. 

15. The method of clarni 10, wherein the first sampling region is substantially at a 
5 center of k-space. 

16. An article of manufacture, comprising: 

a computer readable medium having computer readable program code for 
reconstructing a magnetic resonance (MR) image, the computer readable program 
1 0 code including instructions for causing a computer to: 

obtain an MR data set firom a phased-array coil consistmg of a plurality of 
individual coil-elements and by sampling fu^t and second sampling regions, wherein 
the first region is sampled with a higher sampling density than the second region; 

reconstruct a plurality of images from a plurality of individual coil-elements; 

IS and 

combine the plurality of images fiom the individual coil-elements into a 
resultant MR image using parallel imaging. 

17. The article of manufacture of claim 16, wherem the first sampling region is 
20 near a center of k-space. 

18. The article of manufacture of claim 16, wherein the computer readable 
program code further comprises instructions for causing the computer to fill missing 
k-space locations with zeros before applying a Fourier-transform. 

17 
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19. The article of manufacture of claim 16, wherein flie computer readable 
program code further comprises instructions for causing the computer to apply a 
Fourier-transform. 

5 

20. An apparatus for reconstructing a magnetic resonance (MR) image having a 
number of image pixels, wherein an MR data set is obtained by sampling first and 
second sampling regions, wherein a first region is sampled with a first sampling 
density that is higjier than a second sampling density of a second region, the apparatus 

10 comprising: 

a phased army coil having a plurality of different coil-elements for obtaining 
the MRI data set; 

an MR scaimer for generating MR images corrupted by aliasing artifacts by 
reconstructing the obtained MR data set; and 
15 a processor configured to combine the generated MR images with parallel 

imaging using an acceleration rate of the second region. 

2 1 . The apparatus of claim 20, wherein the MR scanner reconstructs a plurality of 
images fi:om a plurality of individual coil-elements, and wherein the processor 

20 combines the plurality of images fi:om the individual coil-elements into a resultant 
image. 

22. The apparatus of claim 20, wherein the processor inverts a number of 
matrices, the number of matrices being equal to a number of image pixels in the 
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resxiltant MR image divided by the acceleration rate of the second region, wherem 
each matrix has a number of elements equal to or less than the nimiber of coil- 
elements in the phased array coil times the acceleration rate of the second region. 

23. The apparatus of claim 20, wherein the first sampling region is near a center of 
k-space, 

24. The apparatus of claim 20, wherein the MR scanner fills missing k-space 
locations with zeros. 

25. The apparatus of claim 20, wherein the processor is further configured to 
apply a Fourier-transform after filling missing k-space locations with zeros. 
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